Mon. Not. R. Astron. Soc. OOP. [TlfTTI C2007') 



Printed 18 February 2009 



(MN r^TfejX style file v2.2) 



Absolute parameters of the eclipsing binary V821 Cas 
from UBVRI light curves and radial velocities/ f 



0^ 
O 
O 

00 



6. gakirh^, C. ibanoglu^ S. Bilir^, E. Sipahi^ 

^ Ege University, Science Faculty, Department of Astronomy and Space Sciences, 35100 Bornova, Izmir, Turkey 
^Istanbul University, Science Faculty, Department of Astronomy and Space Sciences, 34119 University, Istanbul, Turkey 



Accepted 2008 Month Day. Received 2008 Month Day; in original form 2008 ??? 



6 



> 

cn 
oo 

o 

m 

(N 
O 
On 
O 



ABSTRACT 

We present UBVRI photometric measurements and spectroscopic observations of 
the double-Uned ecUpsing binary V821 Cas. The radial velocities were obtained by 
means of the cross-correlation technique. Simultaneous analyses of the multi-band 
light curves and RVs give the absolute parameters for the stars as: Mi=2.05±0.07 
Mq, M2=1.63±0.06 Mq, Ri=2.31±0.03 Rq, R2 = 1.39±0.02 Rq, Te//i =9 400±400 K, 
and Tef =8 600±400 K. Analysis of the 0-C residuals yielded an apsidal motion in the 

binary at a rate of tj^O .0149±0 .0023 cycle , corresponding to an apsidal period of 
U=118±19 yr. Subtracting the relativistic contribution we find that log fc2o6s=-2.590 
which is in agreement with the value predicted by theoretical models. Comparison 
with current stellar evolution models gives an age of 5.6 x 10^ yr for the system. 



Key words: binaries: stars: close - binaries: eclipsing-binaries: general 
spectroscopic - stars: individual: V821Cas 



binaries: 



1 INTRODUCTION 

V821 Cas (BD +52° 3571, Tycho 4001-1445-1, ¥=8™. 31, (B- 
V)=0.11) was discovered to be an eclipsing binary of early 
A spectral type, with a variability period of ~1.8 days by 
the Hipparcos satellite (ESA, 1997). Degirmenci et al. (2003, 
2007) determined an improved light curve and listed 25 pho- 
tographic times of minima from which they derived an im- 
proved eclipse ephemeris. They obtained BVR light curves, 
and established that the orbit is eccentric from the displace- 
ment of the secondary minimum. Its components do not 
show intrinsic variability, they are well-separated, and well 
within their Roche lobes, which values them proper tests for 
stellar models. 

Several times of minima for V821 Cas have appeared 
in the literature since, but no radial velocity measurements 
and multi-band light curves have been published. The aim 
of this work is to obtain the absolute parameters of the sys- 
tem. Along the next section we first present our spectro- 
scopic data and the set of (75 FT?/ light curves. Making use 
of the spectra we determine the parameters of the radial ve- 



* Based on observations collected at Catania Astrophysical Ob- 
servatory (Italy) and TUBITAK National Observatory (Antalya, 
Turkey). 

I Table 1 is only available in electronic form at the CDS via 
anonymous ftp to http://www.blackwell-syngery.com/doi ... 
I E-mail: omur.cakirli@ege.edu.tr 



locity curve and perform a reliable estimation of the spectral 
types for both components. Then, using the results from the 
spectroscopy together with the our photometric data we an- 
alyze the light curves to obtain reliable solutions that allow 
to determine the rest of the parameters of the system. 



2 OBSERVATIONS 

2.1 Photometric observations 

We report here new photometry of V821 Cas in the Bessell 
UBVRI bands. The photometric accuracy and the phase 
coverage (over 350 observations) are sufficient to guarantee 
a reliable determination of the light curve parameters thus 
permitting a critical evaluation of stellar models. The obser- 
vations were carried out with the 0.40-m telescope located on 
Mt. Bakirhtepe in September of 2007 at the TUBItAK Na- 
tional Observatory (TUG, located in south of Turkey). The 
telescope is equipped with a Apogee Ikxlk CCD (binned 
2x2) and standard Bessel f/SVT?/ filters. 

The instrument with attached camera provide a field- 
of-view of 11'. 3x11'. 3. By placing V821 Cas very close to 
center of the CCD to get highest accuracy, we managed to 
strategically locate the binary on the chip together with two 
other stars of similar apparent magnitude. As comparison 
star we selected BD -1-52° 3575. The check star chosen to be 
BD +52° 3580. Both stars passed respective tests for intrin- 
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sic photometric variability and proved to be stable during 
time span of our observations. 

We collected a total of 1600 points in the UBVRI bands. 
The resulting V-band magnitude differences (var-comp) are 
listed in Table 1 (available in electronic form at the CDS). 
A typical precision of the differential magnitudes is about 
0.008 mag per measurement. Standard IRAfQ tasks were 
used to remove the electronic bias and to perform the flat- 
fielding corrections. The IRAF task imalign was used to 
remove the differences in the pixel locations of the stellar 
images and to place all the CCD images on the same relative 
coordinate systems. The data were analyzed using another 
IRAF task phot with no differential extinction effects taken 
into account given the relative small separation between the 
target and the comparison and check stars. The phase-folded 
light curves for the whole observations are shown in Fig. 5 
with all bands. 



2.2 Spectroscopic observations 

Spectroscopic observations have been performed with the 
echelle spectrograph (FRESCO) at the 91-cm telescope of 
Catania Astrophysical Observatory. The spectrograph is fed 
by the telescope through an optical fibre {UV-NIR, 100 
yitm core diameter) and is located, in a stable position, in 
the room below the dome level. Spectra were recorded on 
a CCD camera equipped with a thinned back-illuminated 
SITe CCD of Ikxlk pixels (size 24x24 ^m). The cross- 
dispersed echelle configuration yields a resolution of about 
20 000, as deduced from the full width at half maximum of 
the lines of the Th-Ar calibration lamp. The spectra cover 
the wavelength range from 4300 to 6650 A, split into 19 or- 
ders. In this spectral region, and in particular in the blue 
portion of the spectrum, there are several lines useful for 
the measure of radial velocity, as well as for spectral classi- 
fication of the stars. 

The data reduction was performed by using the echelle 
task of IRAF package following the standard steps: back- 
ground subtraction, division by a flat field spectrum given by 
a halogen lamp, wavelength calibration using the emission 
lines of a Th-Ar lamp, and normalization to the continuum 
through a polynomial flt. 

Sixteen spectra of V821 Cas were collected during the 
20 observing nights between August 2 and 23, 2006. Typ- 
ical exposure times for the V821 Cas spectroscopic obser- 
vations were between 2400 and 2600 s. The signal-to-noise 
ratio (S/N) achieved was between 70 and 115, depending on 
atmospheric condition, a Lyr (AOV), 59 Her (A3IV), t Psc 
(F7V), HD 27962 (A2IV), and r Her were observed during 
each run as radial velocity and/or rotational velocity tem- 
plates. The average S/N at continuum in the spectral region 
of interest was 150-200 for the standard stars. 



^ IRAF is distributed by the National Optical Observatory, which 
is operated by the Association of the Universities for Research in 
Astronomy, inc. (AURA) under cooperative agreement with the 
National Science Foundation 
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Figure 1. Sample of Cross Correlation Functions (CCFs) be- 
tween V821 Cas and the radial velocity template spectrum (Vega) 
around the first and second quadrature. 



3 SPECTROSCOPIC ANALYSIS 

Double-lined spectroscopic binaries reveal two peaks, dis- 
placing back and forth, in the cross-correlation function 
(CCF) between variable and the radial velocity template 
spectrum as seen in Fig. 1. The location of the peaks allows 
to measure of the radial velocity of each component at the 
time of observation. The cross-correlation technique applied 
to digitized spectra is now one of the standard tools for the 
measurement of radial velocities in close binary systems. 

The radial velocities of V821 Cas were obtained by 
cross-correlating of echelle orders of V821 Cas spectra with 
the spectra of the bright radial velocity standard stars a Lyr 
(AOV), 59 Her (A3IV) and t Psc (F7V) (Nordstrom et al., 
2004). For this purpose the IRAF task fxcor was used. 

Fig. 1 shows examples of CCFs of V821 Cas near the 
first and second quadrature. The two non-blended peaks 
correspond to each component of V821 Cas. We applied the 
cross-correlation technique to five wavelength regions with 
well-defined absorption lines of the primary and secondary 
components. These regions include the following lines: Si ill 
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Figure 2. Radial velocity curve folded on a period of 1.7698 days, where phases 0,1,... are defined to be at mid-primary eclipse. The 
curves show the fit for an eccentric orbit, to the primary and secondary radial velocities. Points with error bars (error bars are masked 
by the symbol size in some cases) show the radial velocity measurements for the components of the system (primary: circles, secondary: 
squares) . 



4568 A, Mgii 4481 A, Hei 5016 A, Hei 4917 A, Hei 5876 
A. The stronger CCF peak corresponds to the more mas- 
sive component that also has a larger contribution to the 
observed spectrum. To better evaluate the centroids of the 
peaks (i.e. the radial velocity difference between the target 
and the template), we adopted two separate Gaussian fits 
for the case of significant peak separation. 

The radial velocity measurements, listed in Table 2 to- 
gether with their standard errors, are weighted means of the 
individual values deduced from each order (see, e.g., Frasca 
et al. 2006). The observational points and their error bars 
are displayed in Fig. 2 as a function of orbital phase as cal- 
culated by means of the ephemeris based on the photometric 
times of the primary eclipse described in Degirmenci et al. 
(2007). 

The first detailed solution of both radial velocity curves 
of V821 Cas components is presented in this study. We found 
the semi-amplitude of the more massive, more luminous 
component to be Ki = 120±2 km s~^ and K2 = 150±2 km 
s"'^ for the secondary component. 



3.1 Spectral classification 

We have used our spectra to classify the primary component 
of V821 Cas. For this purpose we have degraded the spectral 
resolution from 20 000 to 3 000, by convolving them with 
a Gaussian kernel of the appropriate width, and we have 



measured the equivalent width (EW) of photospheric ab- 
sorption lines useful for the spectral classification. We have 
followed the procedures of Hernandez et al. (2004), choos- 
ing hydrogen and helium lines in the blue-wavelength re- 
gion, where the contribution of the secondary component 
to the observed spectrum is negligible. From several spectra 
we measured EWii-y = 11.3 ± 0.9 A, EWhc, = 9.89 ± 0.11 A, 
EWhp = 10.79 ± 0.13 A, and £WMgiiA448i = 0.31 ± 0.04 A. 

From the calibration relations iSW-Spectral-Type of 
Hernandez et al. (2004), we have derived a spectral type AlV 
with an uncertainty of about 1.0 spectral subclass. The effec- 
tive temperature deduced from the calibrations of Drilling 
& Landolt (2000) or de Jager & Nieuwenhuijzen (1987) is 
about 9 400K (AlV). The spectral-type uncertainty leads to 
a temperature error of ATcfr ~ 400 K. 



3.2 Reddening 

The measurement of reddening is a key step in determining 
the absolute temperature scale (and therefore the distance) 
of eclipsing binaries. In addition to moderate distance deter- 
mined by the Hipparcos mission, some reddening is expected 

o 

for V821 Cas due to its low galactic latitude (/=115 .10, 
b= -8° .40). 

Our spectra cover the interstellar Nal (5890 and 5896 A) 
doublets, which is excellent estimators of the reddening as 
demonstrated by Munari & Zwitter (1997). They calibrated 
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Figure 3. Observed spectrum of V821 Gas (large dots) in tlie A4861 region (upper panel) and Mgll A4481 (bottom). The synthetic 
spectrum (A1V+A4V) is displayed with continuous line in the same boxes. The differences (observed-synthetic, shifted) are plotted in 
the bottom of each panel. 



a tight relation linking the Na I D2 (5890 A) and Kl (7699 A) 
equivalent widths with the E(B-V) reddening. On spectra 
obtained at quadratures, lines from both components are 
un-blended with the interstellar ones, which can therefore 
be accurately measured. We derive an equivalent width of 
0.34±0.02 A for only Nal, which corresponds to E(B-V) = 
0.147±0.011 mag. Kl interstellar line is out of our spectral 
range as given in wavelength region in previous section. 



3.3 Rotational velocity 

The width of the cross-correlation profile is a good tool for 
the measurement of i;sin i (see, e.g., Queloz et al. 1998). The 
rotational velocities {v sin i) of the two components were ob- 
tained by measuring the FWHM of the CCF peaks in nine 
high-S/N spectra of V821 Gas acquired close to the quadra- 
tures, where the spectral lines have the largest Doppler- 
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Table 2. Radial velocities of the V821 Cas' components. The 
columns give the heliocentric Julian date, the orbital phase (ac- 
cording to the ephemeris given by Degirmenci et al. 2006), the 
radial velocities of the two components with the corresponding 
errors, and the average S/N of the spectrum. 



HJD 


Phase 


Star ] 


L 


Star 


2 


<5/Ar> 


2 453 000+ 




Vp 


cr 


V, 


(7 




53953.6017 


0.3363 


-93.1 


6.1 


124.8 


9.3 


88° 


53955.6040 


0.4677 


-3.3 


3.2 






76 


53957.5655 


0.5761 


75.0 


6.4 


-88.0 
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80 
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61.0 
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12.1 


75 
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84.0 
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-104.0 


8.1 


95° 


53975.5607 


0.7443 


107.5 


4.8 


-135.9 


8.5 


106° 


53980.5290 


0.5516 


50.6 


4.9 


-73.1 


9.9 


90 


53981.5213 


0.1123 


-75.3 


3.6 


100.8 


7.1 


90° 


53982.5825 


0.7119 


108.0 


2.7 


-127.0 


7.3 


110° 


53983.5405 


0.2533 


-133.0 


2.9 


159.4 


6.6 


113° 


53983.6089 


0.2919 


-121.5 


2.6 


152.4 


3.7 


95 


53984.5390 


0.8175 


100.4 


3.0 


-123.2 


5.2 


97 


53984.6180 


0.8621 


89.3 


1.8 


-104.9 


7.1 


101° 


53985.5492 


0.3883 


-53.4 


3.4 


76.1 


6.0 


100° 


53985.6109 


0.4231 


-40.2 


3.9 


42.8 


5.6 


72 



° Used also for rotational velocities (usini) measurements. 



shifts. In order to construct a calibration curve FWHM- 

vsmi, we have used an average spectrum of HD 27962, ac- 
quired with the same instrumentation. Since the rotational 
velocity of HD 27962 is very low but not zero (wsini ~11 
km s~^, e.g., Royer, Zorec & Fremat 2004 and references 
therein), it could be considered as a useful template for 
A-typc stars rotating faster than t;sini ~ 10 km s~^. The 
spectrum of HD 27962 was synthetically broadened by con- 
volution with rotational profiles of increasing v sin i in steps 
of 5 km s~^ and the cross-correlation with the original one 
was performed at each step. The FWHM of the CCF peak 
was measured and the FWHM-u sin i calibration was estab- 
lished. The vsini values of the two components of V821 Cas 
were derived from the FWHM of their CCF peak and the 
aforementioned calibration relations, for a few wavelength 
regions and for the best spectra. This gave values of 70±1 
km s^^ for the primary star and 57±1 km s~^ for the sec- 
ondary star. 

We classified the spectral type of the primary compo- 
nent as an Al type main-sequence star. The secondary com- 
ponent appears to be an A4V star (see § 5). The spectral 
types of the standard stars a Lyr and 59 Her are very close to 
the primary and secondary component of the V821 Ceis, re- 
spectively. For the construction of the synthetic spectrum of 
the system the spectra of the standard stars, obtained with 
the same instrumentation, have been rotationally broadened 
by convolution with the appropriate rotational profile and 
then co-added, properly weighted by using physical param- 
eters (Ti, T2, 7?i, 7?2, vsini\2) of the components as input 
parameters and, then, Dopplcr-shifted according to the ra- 
dial velocities of the components. In Fig. 3 we show the 
simulation of the spectrum of V821 Cas with the standard 
stars. 



Table 4. Apsidal motion parameters for V821 Cas. 



Anomalistic Period Pa 


1.769813±0.000040 


Reference minimum time To HJD 


51 767.3481±0.0008 


Periastron longitude at Ta,u)° 


148±4 


Apsidal motion rate u) cycle" 


0.0149±0.0023 


Orbital eccentricity e 


0.138±0.011 


Apsidal motion rate U yr 


117±19 



4 EPHEMERIS AND APSIDAL MOTION 

V821 Cas has a moderate eccentric orbit. All available times 
of minimum light covering about 3100 orbital cycles were 
collected by Degirmenci et al. (2007). Their analysis of the 
O-C residuals yields for the first time a period of apsidal 
motion to be 158 yr. 

All available times of minimum light were collected from 
literature are listed in Table 4 together with the new timings 
obtained in this study. The O-C residuals indicate the differ- 
ences between the observed and calculated times obtained 
with linear ephemeris, 

Minl(HJD) = 2 451 767.3504 + 1''.7697397 x E (1) 

The method of Gimenez & Garcia-Pclayo (1983) was 
adopted to carry out the apsidal motion parameters with the 
use of differential-correction approximations. The method is 
successively relates the derived apsidal motion parameters 
to the observed times of minima. We used some subrou- 
tines in the analysis of the O-C values, supplied by O. L. 
Degirmenci. Two secondary and one primary eclipse timings 
given by Otero (2005) were omitted from the analysis due 
to the large deviations from the others. These timings were 
estimated from very few observations, either in ascending 
or descending branch of the eclipses. Assuming the orbital 
inclination i ~ 83° from the light curve analysis (see next 
section) and using the residuals in Table 4 we obtained the 
parameters of apsidal motion with a weighted least-square 
solution. 

The apsidal motion resulting from the fit is a;=0°.0149 
±0°.0023 cycle"^, which is significant at the 6.5a" level. The 
sidereal and anomalistic periods {Ps and Pa), as well as the 
apsidal motion period U — 118±19 yr arc listed in Table 4. A 
plot of the O-C deviations of the times of minimum from the 
linear terms of the apsidal motion is shown in Fig. 4 along 
with the predicted deviations. In the bottom panel of Fig. 
4 we show the differences between observed and computed 
timings, taking into account apsidal motion. The apsidal 
motion parameters given in Table 4 are only preliminary 
results because the data we used cover only 14 % of the 
apsidal motion period. 



5 SIMULTANEOUS LIGHT CURVES AND 
RADIAL VELOCITIES ANALYSES 

In order to obtain the values of the main physical parameters 
of the binary components and to determine the orbital incli- 
nation (i) of the system (necessary to deduce stellar masses 
from the Mi, 2 sin"^ i derived from the RV curves) we applied 
a numerical eclipsing binary model to the UB VRI ohservar 
tions. We choose the Wilson-Devinney (W-D) code imple- 
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Figure 4. Representation of tlie best-fitting apsidal motion parameters. The upper panel shows the observed times of primary (open 
circles) and secondary (filled circles) minima, continuous line are best-fitting curves of primary and secondary minima. The lower panel 
shows the residual of the O-C values from the fitted curves. 



mented into the PHOEBE package tool by Prsa & Zwitter 
(2005) for the LC and differential correction (DC) fits. 

The code was set in Mode-2 for detached binaries with 
no constraints on the potentials. The effective temperature 
(Ti) of the primary star was fixed at 9400±400 K defined in 
§3. The simplest considerations were applied for the param- 
eters of the stars in the model, that is, stars were considered 
as black bodies, approximate reflection model (MREF=1) 
was adopted, and no third hght or spots were included. 
Gravity-darkening exponents gi — g2 — 1 and bolometric 
albedos Albi = Alb2 = 1 were set for radiative envelopes. 
We used square root limb-darkening law. Bolometric limb- 
darkening coefficients were taken from Van Hamme (1993). 

The mass ratio, q=M2/Mi, is very important parame- 
ter in the hght curve analysis, because the WD code is based 
on Roche geometry which is sensitive to this quantity. The 
meiss ratio of 0.795 determined from the radial velocities was 
kept as a fix value. The iterations were carried out automat- 
ically until convergence and a solution was defined as the 
set of parameters for which the differential corrections were 
smaller than the probable errors. The hght curves were ana- 
lyzed individually and the weighted means of the parameters 
i, T2, Qi, fl2, ri and r2 were computed. Our final results are 



listed in Table 5 and the computed hght curves are shown 
as continuous lines in Fig. 5. The uncertainties assigned to 
the adjusted parameters are the internal errors provided di- 
rectly by the Wilson-Devinney code. In the bottom panel of 
Fig. 5 the residuals between observed and computed fluxes 
are also plotted. The residuals reveal that the binary model 
may represent the observed hght curves successfully. 



6 DISCUSSION AND CONCLUSIONS 

6.1 Absolute dimensions and distance to the 

system 

Combination of the parameters obtained from light curves 
and RVs yield the absolute dimensions of the system, which 
are presented in Table 6. The standard deviations of the 
parameters have been determined by 

JKTABSDI^fl 

code, 

which calculates distance and other physical parameters 
using several different sources for bolometric corrections 
(Southworth et al. 2005). The radii of the components were 

^ This can be obtained from 

|http://www.astro.keele.ac.uk/~jkt/codes.html| 
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Table 3. Times of minimum light of V821 Cas. The O-C values refer to the difference between the observed and calculated values. 



Cycle number 


Minimum time 
(HJD-240 0000) 


O-C 


Ref. 


Cycle number 


Minimum time 
(HJD-240 0000) 


O-C 


Ref. 


-2002 


48224 3fiSn 


0.0377 


\ 


672.5 


52957 421 5 


-0.0762 


4 


-1927.0 


OQl n 


0.0303 
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1051.0 




0.0727 
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-1923.5 


4K'^fi'^ 9940 


u. uouo 




1072.5 


O O U U O . O -L 1 U 


-0.0787 


4 


-1846.0 


'-tooyju .'^'-toc! 


0.0363 


2 


1006.5 


0004:0. 0-L±U 


-0.0796 




-lOl -L.O 


"4:0001 .4:10U 




]^ 


1042 


OOUl -L .'±00^ 


n 0791 


7 


-1689.5 


'iO ill. 0^'-t\J 


u. uiyo 


\ 


1262.5 


OIUU J- . ouo^ 


-0.0786 


Y 


-1541.0 


4Qn4.n 91 QO 


u.uooy 


\ 


1272.0 


•^401 S '^'^04 

O^U-LO.OOU^ 


0.0741 


7 


-197.5 


^1417 flQ4n** 


-0.1306 


\ 


1294.5 


^40 ^^i 9n^n 


-0.0704 


g 


-179.0 






\ 


1299.0 


o^uuu.ouy 1 


0.0704 


g 


-162.5 


^147Q fi'^Td** 


-0.1285 


\ 


1299.0 


o^uuu. J- xo 


0.0725 




-26.5 


51720.3851 


-0.0649 


3 


1299.0 


54066.3140 


0.0747 


8 


-26.0 


51721.3976 


0.0627 


3 


1315.5 


54095.3595 


-0.0804 


8 


0.0 


51767.4100 


0.0619 


4 


1315.5 


54095.3601 


-0.0798 


8 


4.0 


51774.4893 


0.0622 


4 


1315.5 


54095.3610 


-0.0789 


8 


17.0 


51797.4962 


0.0625 


4 


1315.5 


54095.3625 


-0.0774 


8 


17.0 


51797.4967 


0.0630 


4 


1315.5 


54095.3636 


-0.0763 


8 


21.5 


51805.3300 


-0.0675 


4 


1471.5 


54371.4442 


-0.0751 


8 


29.5 


51819.4840 


-0.0714 


4 


1471.5 


54371.4452 


-0.0741 


8 


38.5 


51835.4153 


-0.0678 


4 


1475.5 


54378.5192 


-0.0790 


8 


2.30.0 


52174.4524 


0.0643 


4 


1476.0 


54379.5584 


0.0753 


8 


469.0 


52597.4220 


0.0662 


5 


1476.5 


54380.2882 


-0.0798 


8 


630.0 


52882.3561 


0.0723 


6 


1477.0 


54381.3278 


0.0750 


8 


672.5 


52957.4215 


-0.0762 


4 


1490.0 


54404.3347 


0.0753 


8 



** Rejected from the fit owing to a large O-C value. 

Ref: (1) Otero (2005), (2) ESA (1997), (3) Bulut & Demircan (2003), (4) Degirmenci et al. (2003, 2007), (5) Ak & Filiz (2003), (6) Baki§ 
et al. (2003), (7) Brat et al. (2007), (8) This study. 



Table 5. Results from the simultaneous solution of UBVRI hand 
Ught curves of V821 Cas. 



Table 6. Fundamental parameters of V821Cas. 



Parameter 


Value 


in 


82.6±0.1 


e(°) 


0.127±0.007 


a;(°) 


155±4 


Ti (K) 


9400 [Fix] 


T2 (K) 


8600±17 


Qi 


5.067±0.014 


O2 


6.693±0.035 


Qspec 


0.795 


Li/{Li+2)u 


0.825±0.001 


L\/(L\+2)b 


0.794±0.003 


Li/{L\+2)v 


0.801±0.004 


Li/{Li+2)r 


0.790±0.004 


Li/{Li^2)i 


0.781±0.008 


n 


0.2434±0.0013 


r2 


0.1466±0.0017 


A<j) 


-0.0424±0.0005 




0.002 



estimated with uncertainties of 1.2 % and 1.6%. However, 
the uncertainties on the masses, being 3.3 % and 3.6 %, 
are slightly larger than the criteria of Andersen (1991) for 
the necessary precision of absolute dimensions of stars to be 
used for comparison with theoretical models. 

An inspection of the temperatures, masses and radii of 
the component stars reveals a binary system composed of 
two main-sequence stars. The temperature T^ffi = 9400 
K, mass Mi = 2.04 M© and radius Ri = 2.31 Rq of the 



Parameter 


Primary 


V821 Cas 

Secondary 


Spectral Type 


AlVitl 


A4V±1 


a (Rq) 




9.496±0.099 


(km s-i) 




-0.05±0.01 


in 




82.6±0.1 


1 




0.795±0.017 


Mass (M0) 


2.046±0.067 


1.626±0.058 


Radius (R0) 


2.311±0.028 


1.392±0.022 


log 9 {cgs) 


4.021±0.007 


4.362±0.012 


T.ff (K) 


9400±400 


8600±400 


{vsin i)obs (kni s ^) 


70±1 


57±1 


{vsin i)calc. (km s~^) 


66itl 


40±1 


log (L/Lq) 


1.58±0.06 


0.98±0.01 


d (pc) 




260±12 


J, H, Ks (mag)* 


7.978±0.018, 


8.007±0.027, 7.956±0.027 


ficCosS, /i^(mas yr"-*-)** 


-2.51±0.59, -5.95±0.62 


U,V,W (km s-l) 


4.77±0.70, 1.07±0.34, -6.82±0.80 



*2MASS All-Sky Point Source Catalogue (Cutri et al. 2003) 
**Newly Reduced Hipparcos Catalogue (van Leeuwen 2007) 



primary are consistent with the spectral type of AlV, and 
the temperature Teffi — 8450 K, mass M2 = 1.62 M0 and 
radius R2 — 1.39 Rq of the secondary are consistent with 
an A4V spectral type star. 

The colour excess E{B — V) for a star may also be 
determined from the wide-band photometry. For calculation 



8 0. (^akirh et al. 




Figure 5. Observed phased light curves of V821Cas with the best-fitting Phoebe model light curves. For clarity the UBVRI light curves 
and residuals of the fit (lower panel) have been offset by constant amounts as shown in figure. The right panel expands the region near 
the primary minimum for the R-band. 



the reddening-free Q — parameter we used the well-known 
equation below, 

Q = {U-B)~ [E{U - B)/E{B - V)] x {B - V), (2) 

for stars from 09 to A2. The average ratio of colour excesses 
for stars from 08 to A2 was adopted as [E{U — B)/E(B — 
V)] = 0.72 ± 0.03 (Hovhannessian 2004). The standard re- 
lation between Qubv and {B — V)o was also given in Hov- 
hannessian's study. For calculation of the total absorption 
in the visual magnitude the following relation, the ratio of 
selective-to-total extinction in the V band, 

[Av/E{B - V)] = 3.30 -I- 0.28(B - V)o + OME{B ~ V), (3) 

given by Drilling & Landolt (2001) was adopted. 

Using the observed visual magnitude and color indexes 



given by Oja (1985) as V=8.26, U-B=0.07 and B-V=0.11 
mag we compute the reddening-free parameter as Q = 
—0.009. Using the tables given by Hovhannessian (2004) we 
find the intrinsic colour index of (B — V)o — —0.021 mag 
and then the colour excess of E{B — V) = 0.131 ±0.020 mag 
for the system. The colours of the system were the aver- 
age of three measurements. We should mention that there is 
no clue about the date of observations obtained, i.e. orbital 
phase. This value is in agreement with the E(B-V)=0.147 
determined from the spectra given in §3.2. 

Using the two E{B — V) values derived from photomet- 
ric and spectroscopic data we calculated the de-reddening 
distance modulus of the system. To estimate the bolomet- 
ric magnitudes of the components we adopted Mi,oi=4.74 
mag for the Sun. Using the bolometric corrections given by 
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Drilling & Landolt (2000) and Girardi et al. (2002) we esti- 
mate the distance to the system as 260 and 274 pc, respec- 
tively, with an uncertainty of 12 pc. However, the average 
distance to the system is estimated to be 2O6I29 P*- from the 
trigonometric parallax measured by the Hipparcos mission. 



6.2 Internal structure 

Our detection of the apsidal motion rate of the V821 Cas in 
§4 provides the opportunity to test the models of stellar in- 
ternal structure. In addition to the classical Newtonian con- 
tribution, the observed rate of rotation of apsides includes 
also the contributions arising from General Relativity. The 
theory of General Relativity estimates the relativistic con- 
tribution to observed rate as the following Einstein formula: 



LOrel — 5.45 X 10 



1 , Ml+M2 2/3 
1 - e2 ^ Pa ' 



(4) 



where Mi(i = 1,2) denotes the individual masses of the com- 
ponents in solar mass, e is the orbital eccentricity (from ra- 
dial velocity analysis) and Pa is the anomalistic period of the 
system in days. Use of equation (4) for V821 Cas yields the 

o 

relativistic contribution of .0009 cycle^^ which is only 6 
per cent of the observed rate. Removing the relativistic con- 
tribution from the observed rate of apsidal motion, one can 
calculate the observational average value of internal struc- 
ture constant {k2obs) using the following formula: 



^2063 = 



C21 + C22 360 ' 



(5) 



where C2i(* = 1,2) are the functions of the orbital eccen- 
tricity, fractional radii and masses of the components and 
ratio of the rotational velocities of the component stars to 
the Keplerian velocity. The observed average value of log k2 
is found to be -2.56±0.07 from Eq. (5). The theoretical in- 
ternal structure constants for the components (log k2theoi,2) 
are taken from the theoretical calculations of Claret (2004). 
Among the tabulated values, interpolation for the masses 
and log g of the component star V821 Cas yielded -2.48 and 
-2.37 for the primary and secondary components, respec- 
tively. The mean theoretical value of internal structure con- 
stant was computed as log fc2t;ieo=-2.43 which is in agree- 
ment with that of observed within 2a level. The internal 
structure constant obtained from the 0-C analysis is smaller 
by 35 % than predicted from the theory. The present anal- 
ysis shows that the component stars are more concentrated 
in mass than predicted by theoretical calculations. 



6.3 Evolutionary stage and age of the system 

In Fig. 6, we plot the location of V821 Cas stellar compo- 
nents in log Te/ / — log g diagram. The evolutionary tracks 
for masses 1.6, 2.0 Mq and z = 0.008 (Girardi et al. 2002) 
are also shown in this figure. Higher metallicity tracks given 
by Girardi et al. do not match with the observed properties 
of the components. An isochrone corresponding to an age 
of 5.6x10* yr, computed from Girardi et al. (2002)0 is also 
shown in Fig. 6. The effective temperature of the primary 



^ We computed the Isochrones using the facility available in the 
web site: http: / /stev.oapd.inaf.it/~Igirardi/cgi-bin/cmd| 



component appears to be slightly higher than the model for 
2 Mq . However, when the uncertainties of the effective tem- 
peratures of the components are taken into account it is seen 
that the deviations from the models are less than la level. 
The diagram suggests that both components of V821 Cas are 
main sequence stars, relatively young with an age of about 
560 Myr. 



6.4 Population type and kinematical analysis of 
V821 Cas 

To study the kinematical properties of V821 Cas, we used 
the system's center-of-mass' velocity, distance and proper 
motion values, which are given in Table 6. The proper mo- 
tion data were taken from newly reduced Hipparcos cata- 
logue (van Leeuwen 2007), whereas the center-of-mass ve- 
locity and distance are obtained in this study. The system's 
space velocity was calculated using Johnson & Soderblom's 
(1987) algorithm. The U, V and W space velocity compo- 
nents and their errors were obtained and given in Table 6. 
To obtain the space velocity precisely the first-order galac- 
tic differential rotation correction was taken into account 
(Mihalas & Binney 1981), and -1.08 and 0.65 kms"^ differ- 
ential corrections were applied to U and V space velocity 
components, respectively. The W velocity is not affected in 
this first-order approximation. As for the LSR correction, 
Mihalas & Binney's (1981) values (9, 12, 7)© kms"^ were 
used and the final space velocity of V821 Cas was obtained 
as S = 19.36 kms~^. This value is in agreement with other 
young stars space velocities. 

To determine the population type of V821 Cas the 
galactic orbit of the system was examined. Using Dinescu 
et al. (1999) N-body code, the system's apogalactic (Rmax) 
and perigalactic (Rmin) distances were obtained as 9.09 and 
8.07 kpc, respectively. Also, the maximum possible vertical 
separation from the galactic plane is |2maa;|=50 pc for the 
system. When determining the ellipticity the following for- 
mula was used: 



Rn 



Rn 



Rn 



+ Rn 



(6) 



The ellipticity was calculated as e = 0.059. This value shows 
that V821 Cas is orbiting the Galaxy in an almost circular 
orbit and that the system belongs to the young thin-disc 
population. 
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Figure 6. Location of the two stellar components of V821 Cas in log T^ff- log g diagram, together with evolutionary models for the 
masses of 1.6 and 2.0 (Girardi et al. 2002) and the isochrones of an age t=5.6xl0* years (long-dashed-dotted line). The open circle 
corresponds to the primary and the open square to the secondary. 
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